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ABSTRACT: 


A f'rst  principles  kinetic  theory  is  developed,  applicable  to  annular  dependent 
emission  in  Electron  Stimulated  Desorption.  Expressions  for  the  ionic  and  neutral 
atom  ESD  cross  sections  arc  formulated  and  applied  in  a model  calculation  o?  O 
<•.  e ion  from  W(ijl),  Stiong  focussing  of  the  outgoing  ions  was  found,  with  the 
use  of  a modei  ion-solid  potential  in  which  the  substrate  was  free  ot  excitation, 
f . spot  groups  were  simulated.  The  peak  ion  energies  obtained  with  this  mode! 

however,  small  compared  to  experimental  energies  for  the  high  coverage  case. 
The  need  to  introduce  substrate  excitations  to  describe  this  case  is  discussed. 
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I.  INTRODUCTION 

Much  work  has  been  done  in  recent  years  on  the  phenomenon  of  Electron 
Stimulated  Desorption  (ESD) . In  these  experiments,  a beam  of  low  energy 
electrons,  in  the  range  of  100-500  ev  impinges  on  an  adsorbate  covered  sur- 
face, causing  the  ejection  of  ions  or  neutral  particles.  Several  articles  re- 

1-4 

viewing  ESD  have  appeared. 

Experimental  data  have  been  collected  about  the  ion  yields  and  cross 

5 

sections.  Energy  distributions  of  emerging  ions  have  been  measured.  In  the 

last  several  years,  the  angular  distribution  of  ejected  ions  has  been  investi- 
6-8 

gated.  Strong  focussing  of  the  desorbing  ions  has  been  revealed.  Desorption 
patterns  were  obtained  which  contain  sets  of  well  defined  cones  for  ion  emis- 
sion, both  normal  and  non-normal  to  the  surface.  The  number  and  orientation 
i of  emission  cones  is  related  to  the  substrate  geometry. 

ESD  Ion  Angular  Distributions  may  thus  provide  valuable  clues  to  the  ad- 
sorption geometry,  on  the  scale  of  the  unit  cell  size.  This  direct  information 

■ 

is  complementary  to  that  obtained  , for  example,  by  LEED,  which  is  sensitive 

I ■ 

| « 

tc  mainly  long  range  order  in  the  adsorbed  layer  and  substrate. 

R ' 

In  this  theoretical  work,  we  address  ourselves  primarily  to  the  angular 

dependence  problem.  Previous  theoretical  treatments  of  ESD  have  been  strict - 

9,10 

ly  one  dimensional,  and  phenomenological  in  flavor.  This  work  presents  a 

I i 

first  principles  transport  formulation  of  the  desorbing  ion  dynamics.  The 
ion  -solid  interactions  are  assumed  to  be  time  Independent,  or  long  lived 
relative  to  ion  flight  times  . 


2 

Section  n establishes  the  conceptual  framework  and  introduces  kinetic 
equations  which  can  be  solved  formally,  in  three  dimensions,  for  the  linear 
response  to  external  electron  bombardment.  In  Section  III,  expressions  for  the 
ionic  and  neutral  atom  ESD  cross  sections  are  obtained  from  the  kinetic  theory  , 
in  a form  suitable  to  efficient  numerical  integration.  Section  IV  is  an  illustra- 
tive model  calculation  of  ESD  Ion  Angular  Distributions . The  experimentally 
interesting  case  of  0+  desorption  from  W (111)  was  chosen,  with  the  surface 
assumed  free  of  excitation  during  the  ion  ejection.  Section  V is  reserved  for 
discussion. 

We  found  that  ion  interactions  with  a ground  state  surface  can  replicate  the 
focussing  and  ion  distributions  in  experimental  ESD  patterns,  for  several  ad- 
sorption geometries.  However,  this  model  surface  potential  probably  cannot  ex- 
plain the  large  0+  ion  energies  observed  at  high  coverage.  We  speculatively 
suggest  that  long  lived  excitations  on  the  surface  may  play  the  doinlnent  role 
in  ESD  of  0+  at  high  coverage . 


,$*v 
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II.  FORMULATION  OF  THE  PROBLEM 

This  calculation  employs  a multi-step  picture  of  the  desorption  process 

1-4,9-11 

which  has  been  found  useful  in  previous  ESD  investigations.  Absorbed  particles 
are  Franck-Condon  excited  to  an  anti-bonding  state  from  which  ionic  desorption 
may  occur.  Propagation  follows,  with  the  dynamics  determined  by  an  effective 
adparticle-surface  interaction  potential.  Recapture  may  occur  via  bond  healing 
transitions.  Reneutralized  particles  may  still  possess  sufficient  kinetic  energy 
to  escape  as  neutrals,  in  which  case  ground  state  propagation  completes  the 
desorption  process. 

The  above  picture  is  useful  for  describing  the  response  of  adparticle  sub- 
strate systems  which  are  interacting  weakly  with  an  external  beam.  Hence,  we 
linearize  the  kinetic  equations  presented  below  in  the  incident  flux.  Adparticle 
propagation  is  described  classically.  The  excitation  and  decay  mechanisms  in- 
volve electronic  transitions,  and  thus  require  quantum  mechanical  descriptions. 

A single  atom  or  molecule  interacting  with  a surface  may  be  described 
by  the  following  Hamiltonian: 

H = -■£  + Hje.v,) 

£r\  (D 

Here,  Hg  contains  all  the  potentials  in  the  system,  as  well  as  kinetic  energy 
operators  for  electrons  which  participate  in  atom-surface  interactions.  The 
kinetic  energy  operator  shown  acts  on  adparticle  mass  center  coordinates,  with 
M denoting  the  adparticle  mass. 

In  the  Born-Oppenheimer  Approximation,  the  propagation  phase  of  de- 
sorption decouples,  to  order  me^M,  from  the  electronic  states  of  the  system. 
Electronic  responses  to  shifts  in  the  adparticle  coordinates  may  be  regarded  as 
instantaneous,  permitting  a set  of  effective  interaction  potentials  to  be  defined. 


Excitation  and  reneutralization  are  likewise  regarded  as  sudden,  inelcstic 
events  in  which  the  internal  adparticle  state,  and  thus  the  propagation  dynamics, 
change  discontinuously , V/e  assume  that  the  asymptotically  recognizable  ad- 
particle  states  partially  diagonalize  H close  to  the  surface. 

For  the  present,  consideration  is  restricted  to  a two  level  system,  with 
n= i for  the  neutrai  ground  state  For  positive  ions,  n=2*  The  z-dependent  be- 
havior of  the  effective  potentials  E^  (?)  is  sketched  qualitatively  in  Fig.  1. 

The  minimum  of  E(f?)  lies  at  the  chemisorption  site.  In  order  for  ionic  desorp- 
tion to  occur,  the  antibond ingcurve  i: ^ (r)  must  have  the  generally  repulsive 
character  shown.  One  expects  E^  t°  become  attractive  at  long  range  due  to 
image  charge  density  residing  on  the  surface.  The  difference  (F 2<f>  - El<3 
is  regarded  as  the  internal  ionization  energy  of  the  adsorbate.  At  large  ad- 
sorbate to  surface  separations  it  is  identical  to  the  free  atomic  value.  Varia- 
tions of  E (?)  in  the  plane  parallel  to  the  surface  are  crucial  for  the  angular 
n 

dependent  desorption  problem,  and  are  discussed  below. 

This  picture  of  desorption,  and  the  transport  theory  developed  below, 
may  be  readily  generalized  to  include  multiple  excited  state  systems,  in  which 
substrate  excitations  and  neutral  excited  states  of  adparticles  must  be  con- 
sidered. In  addition,  photodesorption  may  be  treated  within  basically  the  same 
12 

formalism . 

Little  attention  in  the  literature  has  been  given  to  the  ionic  states  which 
govern  the  desorption  kinetics.  Most  work  has  focussed  on  the  ground  state 
potential  minima,  where  adsorption  occurs.  Studies  of  electron  dissociation  of 
molecules  furnish  a guide  to  the  antibonding  states,  however,  the  model  adopted 
for  in  this  preliminary  study  is  likely  to  be  crude.  The  difficulties  involved 
in  modelling  E-,  for  positively  ionized  cores,  are  less  formidable  than  for  neutrals 
or  negative  ions,  for  which  exchange  and  correlation  effects  may  be  Important. 
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The  functions  En(r)  are  assumed  known  in  developing  a kinetic  theory  below. 
Furthermore,  the  solutions  of  the  classical  equations  of  motion  for  ions  and 
neutral  particles  are  assumed  to  be  known,  and  to  adequately  describe  propaga- 
tior>  oi  the  massive  adparticles.  The  classical  solutions'r*n(t)  depend  para- 
metrically on  time  and  on  the  initial  coordinates  for  each  adparticle  trajectory. 

This  classical  description  is  most  applicable  to  unbound,  and  hence  de- 
sorbing states.  Ions  in  state  2 quickly  acquire  substantial  kinetic  energy. 

The  same  argument  applies  to  atoms  desorbing  in  the  ground  state  after  suffer- 
ing reneutralization.  A classical  description  is  also  acceptable  for  vibrational 
states  of  massive  adatoms  at  high  temperatures,  that  is,  when  ~f>  « kT 

Here,  h<u0  may  be  taken  to  be  the  level  spacing  for  a harmonic  oscillator 
approximation  to  near  the  chemisorption  site. 

When  ~p)  cJ0  £ kt  , classical  kinetic  equations  still  provide 
qualitatively  correct  results  for  the  linear  response.  The  propagation  of  partieles 

only  in  unbound  trajectories  need  be  followed.  The  principal  correction  intro- 
i 

duced  by  quantizing  the  vibrational  states  is  to  modify  the  shape  of  ion  energy 
distributions . 

We  denote  by  R(r),  Q(r)  the  transition  rates  for  ionization  and  reneutrali- 
zation respectively.  These  functions,  assumed  known,  of  adatom-coordinate 

» 

r^are  averaged  over  all  initial  and  final  states  of  the  substrate  and  incident 

’ 

electron  consistent  with  a change  of  adatom  state.  Inasmuch  as  the  dominant 

i 

impulse  to  adsorbed  or  desorbing  particles  is  provided  in  this  picture  by  the 
gradients  of  E the  small  recoil  suffered  during  ionization  or  non-radiative 
decay  will  be  neglected. 

Electronic  energies  generally  employed  in  ESD  experiments  lie  in  the 
range  100  - 500  ev.  So  long  as  incident  electrons  are  of  high  but  non-re lativistic 


6 


k 


energy,  the  excitation  phase  of  ESD  may  be  described  in  tne  Born  Approximation. 
For  simple  ionization  processes,  the  coordinates  of  two  final  state  electrons  are 
averaged  over  in  K(n.  These  exit  quickly  from  the  surface  region,  and  are  either 
absorbed  by  the  solid  or  emitted  as  Auger  electrons  to  the  vacuum.  Inasmuch  as 
incident  electron  energies  are  large  compared  to  (E^  - E^) , ionization  can  occur 
from  within  the  entnc  f ranck- ..  ondcn  region. 

Non-radiative  decay  channels  dominate  transitions  back  to  the  ground  or 
iow  lying  excited  states.  Optical  decays  are  discounted,  as  their  lifetimes  are 
long  compared  to  typical  ionic  flight  times.  Ion  reneutralization  involves  matrix 
elements  for  electron  capture  irom  the  solid.  Hence,  this  process  is  efficient 
at  close  range,  and  weakens  rapidly  with  increasing  distance  from  the  surface. 

The  surface  plane  is  assumed  to  be  tessellated  by  a periodic  grid  of  sur- 
face unit  cells,  each  of  whi'  ' anisotropic  in  the  x-y  direction.  We  regard 
this  structure  as  known  ig  of  the  surface,  or  rearrangements  of  the 


substrate  layer  spacin'  are  known  to  occur  under  certain  conditions, 

can  be  incorporated  into  later  calculations  as  required. 


The  symmetry  and  periodicity  of  the  surface  mesh  are  imposed  on  the 

functions  E , R,  Q which  describe  the  stages  of  desorption.  Transverse 
n 

anisotropies  in  the  adparticle-surface  interactions  become  neglible  while  the 
atom  to  surface  separation  is  small  compared  to  the  sample  dimensions. 

We  formulate  a classical  transport  theory  entirely  within  a single  surface 
unit  cell,  extended  in  the  z direction  normal  to  the  surface.  At  macroscopic 
distances,  each  such  cell  behaves  as  a point  source  of  angular  dependent 
intensity.  Transport  equations  linear  in  the  excitation  rate  R(f)  are  now  obtain- 
ed. The  rate  Q fr)  is  treated  as  a zero  order  quantity.  Later  calculations  verify 
that  the  desorption  cross  sections  depend  exponentially  on  0 (r) . 
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Let  r i d v be  the  piobability  of  finding  an  atom  or  ion  in  state 

n at  time  t,  within  the  phase  space  volume  element  atr  v',  where  these  uie 
the  i^nic  position  and  velocity.  The  following  normalization  is  assume!. 


for  the  present: 


( d ^ dv  fc(\r 


v,t)  = 1 


r - f 2 * 

r>  o-a nr 


This  condition  is  relaxed  after  linearizing  the  theory. 

All  spontaneous  desorption  processes  are  neglected  for  the  systems  we 
consider  at  the  exper >m entail / interesting  temperatures.  Thermal  desorption 
may  occur  as  a by-product  of  local  surface  heating  by  the  incident  electron 
beam.  This  effect  may  be  disregarded  for  weak  fluxes. 

In  the  absence  of  external  perturbation,  the  atom-surface  system  is  taken 
to  be  a gas  of  non-interacting  particles  in  thermal  equilibrium  with  the  solid. 

State  2 is  unpopulated,  hence: 

/f  ( r~y)  ~ 0 

Henceforth,  zero  superscripts  denote  equilibrium  state  quantities.  The  Boltzmann 
equation  which  describes  the  ground  state  is; 

v-tf-  &%-£,•%  f =MJ  - o 

[dtJ5uRf  ( 2 ) 

Here,  M is  the  adsorbate  mass.  The  surface  interaction  term  describes  collisions 
between  adsorbate  atoms  and  substrate  excitations.  In  the  present  case,  these 
provide  principally  a relaxation  mechanism  for  the  ground  state. 

Now  we  allow  the  system  to  be  driven  off  equilibrium  by  an  external 
electron  beam.  Competition  between  the  relevent  transition  mechanisms  is 
summarized  in  the  following  kinetic  equations: 

ff  ~(§l 


f 


Here,  the  subscripts  Q,  E refer  to  the  quenching,  excitation  mechanisms.  We 
made  the  choice  (dP,/dt)  ,.**  0 , in  writing  the  above,  thus  neglecting 
scattering  and  energy  loss  of  the  desorbing  ions  due  to  interaction  with  sub- 
strate excitations.  Such  processes  will  be  studied  in  future  investigations. 

Ions  entering  state  2 are  described  by  the  source  term  (dP^/  dt)j,.  Inas- 
much as  recoil  is  neglected,  one  may  approximate: 


fiH)  ^ R(r,z.-t)R(t-) 

Id-tJe 


The  corresponding  ground  state  depletion  occuring  via  the  ionization  pro- 
cess given  by: 

(4S.)  = - ~ -f4S 


'd-t/a 


Depletion  of  the  ionic  state  through  reneutralization  is  given  by: 

= - fi(W)Q<?) 

( dtJo. 


Reneutralized  adparticles,  which  re-enter  the  ground  state  are  the  only 
ones  which  can  desorb  in  state  1.  These  contribute  a distinctly  non-thermal 
component  to  Pj,  which  must  be  proportional  to  the  following  re  population  term, 
neglecting  the  recoil  upon  quenching: 

g<n 


The  following  set  of  kinetic  equations,  linear  in  R,  describes  the  angular 


dependent  ESD  problem  for  the  weakly  perturbed  case: 

d t 


(3a) 


f?'  = HQ 

d-t 


,1 


(3b) 


All  terms  here  are  of  order  R , as  required.  One  dimensional  analogs  of  Eqns. 

(3a,  3B  were  employed  in  previous  work  on  photodesorption. 

o 

In  obtaining  the  result  above,  we  recalled  that  P^jeO.  Hence,  to  lowest 
order  P-,  P^**  R.  The  ground  state  distribution  was  taken  to  be; 

V'?  ♦ V ’ ' 


D 


where  the  depleted  equilibrium  component  is  defined  to  satisfy: 

■D  P ^ / / o \ p 

l?ut?f  / 


/f  V • % v°-  fifi) 

d-t  i 


The  remaining  contribution,  P^,  is  then  the  solution  to  Eqn.fib).  Inasmuch 
_ D 

as  Pj  satisfies  an  equation  independent  of  Q,  a complete  description  of  ground 

1 n 0 

state  desorption  is  yielded  by  P^  . Eq.(2)was  used  to  approximate  P^  Pj 
to  dominent  order  in  R.  For  weak  coupling  to  the  electron  flux,  the  characteristic 
survival  time  in  the  ground  state  against  ionization  is  long  compared  to  the 
time  for  substrate  induced  relaxation.  The  following  condition  then  holds: 


>7  (dR  ) 
l d*/ £ 


i 


Depletion  of  P®  may  then  be  neglected,  as  its  effect  on  P^1  and  Pj  is  of 
2 

order  R at  least. 
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In  experimental  situations,  one  ideally  observes  the  current  of  ions  or 
neutral  atoms  ejected  in  each  state.  Linearity  of  ESD  ion  currents  with  the 
incident  electron  flux  has  been  supported  by  observation.  Hence,  the  inter- 
mediate and  strong  coupling  limits  of  this  theory  appear  to  be  of  little  prac- 
tical interest  at  present. 

III.  DESORPTION  CROSS  SECTIONS 

The  total  cross  section  for  desorption  in  the  n th  state  is  given  by: 

^ / F&ct) 

Here,  I (z,t)  is  the  net  adparticle  current  in  state  n crossing  a plane  parallel 
to  the  surface  at  z,  at  time  t,  per  unit  cell.  ^ is  evaluated  for  L large  com- 
pared to  unit  cell  dimensions.  In  is  proportional  to  both  the  electron  flux  F, 
and  to  e>(f)  , the  number  of  adsorbed  atoms  per  unit  cell  remaining  on  the 

surface  at  t.  The  cross  sections  are  independent  of  F for  linear  response.  The 
coverage  (Jj)  (fi  and  current  _Z-n  ( L-t  ~t.)  have  the  same  time  dependence 
via  ^ ; hence,  is  time  independent  as  well.  Here  is  the 

desorption  cross  section  summed  over  all  channels,  while  t is  the  elapsed  time. 

For  angular  dependent  ESD  studies,  the  differential  desorption  currents 
1^  (A,  (o  , t) , and  the  differential  cross  sections  are  of  interest. 

Here  Q,  specify  the  asymptotic  emission  directions.  The  asymptotic  speeds 
are  integrated  over.  The  total  and  differential  cross  sections  are  related  by: 

<r„  = 


One  may  find  explicit  expressions  for  <TnC(9}f))  in  terms  of  "/Re- 
distribution functions  by  straightforwardly  following  the  microscopic  adparticle 
flux  within  one  unit  cell: 
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Here  lies  in  the  plane  parallel  to  the  surface.  Aq  is  the  unit  cell  area. 

One  first  evaluates  the  current  density  contribution  frum  d v; 

f,-c)<<3v  « & W l -y)  vd  V 

associated  v/ith  the  phase  point*?,  v. 

In  view  of  the  complexities  associated  with  the  angular  dependence 
problem,  we  adopt  an  approach  well  suited  to  numerical  simulation  of  ESD 
Ion  Angular  Distributions.  The  explicit  result  above  for  the  differential  cross 
sections  is  cumbersome  to  implement  efficiently  on  the  computer. 

The  steady  state  ESD  response  Is  found  by  solving  an  equivalent  time 
dependent  problem.  We  Simulate  the  flux  of  Ions  or  neutral  atoms  desorbed  in 
response  to  an  incident  electron  pulse  of  short  duration  r.  then  take  the  limit 
r -»0  • Suitable  expressions  for  the  total  desorption  cross  sections  are  con- 
structed using  the  transport  formulation.  Numerical  differentiation  in  the 
asymptotic  plane  allows  one  to  recover  the  angular  dependence. 

This  strategy  is  implemented  using  the  time  ordered  picture  of  ESD.  An 
ensemble  of  adsorbed  ground  state  atoms  is  perturbed  at  t * 0 by  the  electron 
beam.  The  forward  evolution  of  the  system  is  followed,  without  further  inter- 
action with  the  external  beam,  to  a time  t,  when  all  ions  or  reneutralized  atoms 
which  can  desorb  have  crossed  a detection  plane.  Each  trajectory  must  be 
individually  followed  to  its  asymptotic  direction  via  the  classical  equations  of 
motion.  One  evaluates  the  contribution  made  by  each  ion  or  atom  to  the  cross 
section.  Including  attenuation  due  to  quenching. 

Time  dependent  solutions  of  the  kinetic  equations  are  required.  The 

formal  solutions  for  P map  the  evolution  of  an  initial  distribution  function  for- 

n 

ward  in  time.  These  are  used  to  construct  the  desorption  cross  sections. 
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By  convention,  the  dependence  of  all  quantities  on  r,  v,  is  suppressed, 
with  only  their  parametric  time  dependence  displayed.  The  classical  dyna- 
mics problem  will  be  solved  numerically.  Liouvilles  Theorem  prompts  the 
replacement  in  Eqns  (3a,  3b): 

_ d 


-f  V oV.  ~t  CL  o 

3* 


Here,  a is  the  acceleration  in  the  n th  state  given  by: 

a = Vf.  £n(r) 

Equations  (3a, 3b)  now  take  the  form: 


d f 1<H  -=  R«.) 

d-t 

d%'  =• 

dt 


(4a) 


(4t) 


First,  the  ionic  survival  probability  against  reneutralization,  S^t^t), 
is  defined  to  be: 


S Jt,t)  = eVL- 

'o 


(5) 


The  exponent  is  the  reneutralization  rate  integrated  along  the  ionic  trajectory 
from"r(t  ) to*r(t),  which  is  governed  by  E2(r).  An  alternative  notation  which 


t,  ’ 
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It  is  useful  to  interpret  S.,  (t  , t)  as  ari  evolution  operator,  which  maps 
ions  from  their  initial  to  final  coordinates.  One  obtains  the  form  of  Fq.(5)ty 
solving  Eq.  (4a)  with  the  ionization  term  neglected. 

ESD  cross  sections  are  experimentally  found  to  be  several  orders  of 
magnitude  smaller  than  atomic  or  molecular  ionization  cross  sections  would 
indicate.  Hence,  quenching  is  presumed  to  be  efficient  and  we  expect 
s2  (t  , t)  <.<1,  when*rttQ)  is  close  to  the  surface,  and"?(t)  lies  in  the  detection 
plane.  It  is  also  clear  that  S7  (t  ,t).S  (t,t‘)  = S,  (t  ,t'),  and 


S,(t  .t)  = S, 

Ca  O “ 


-l 


Ur  t y 


(6a) 


The  following  trial  solution  for  P^1  (t)  is  adopted: 

p21  (t)  : g(t)  s2(to,t) 

After  formally  differentiating  the  above  and  comparing  with  Eq.  (4a)  the  unknown 
function  g(t)  is  found  to  satisfy: 

dg  it')  -=  R(rJ  f^’ct)  (t*,  t') 

The  formal  solution  P2*  (t)  for  the  ionic  state  is 

- g'wsjtit) 

fo 

Here,  we  integrated  the  equation  for  g(t)  formally  over  tQ  <t'  <t.  The  boundary 
condition  g(t  ) * P?*  (t  ) is  found  by  taking  the  limit  t -*t . + €.  One  notes  that 

Lim  S,(t  ,t)  = 1.  Both  R(t) , Q(t)  are  assumed  to  be  mathematically  well  behaved. 

2 ° o o 

The  equilibrium  distribution  satisfied  P^  (t)  * P^  (tQ) . 

The  solution  P^ 1 (t)  yields  the  ionic  distribution  function  which  has 

evolved  from  the  initial  distributions  P2l(tQ)  and  P °(t  ).  The  first  term  in 

eq.  (6a)  attenuates  the  initial  ion  density,  while  the  second  incorporates 

ionizations  which  occur  during  the  Interval  Jt^, t|. 


> V .?VSR?&  V.'.al  i..  ,c.  *• 


The  formal  solution  of  Eq.  (4a),  containing  the  ground  state  component 


of  quenched  particles,  is: 


(6b) 


Here,  one  uses  the  state  l dynamics  to  map  r (t),  v (t* ) forward  to  time  t in  the 
second  term. 

We  construct  the  cross  sections  after  first  examining  the  response  of 
a single  ground  state  adatom  to  ionizing  transitions  during  the  short  finite 
interval/"—  / / & J,  Let  U+(-r ,0)  be  the  probability  of  finding  the  atom  as 

an  ion  at  t=0.  Using  Eq.  (6a)  this  is  found  to  be; 

U*('T,o)  - f 'eft' ft (fJ SJt'.o) 

/•r  (7) 


Here,  the  choice  P2*  (-r  ) * 0 was  made,  inasmuch  as  the  initial  ionic  population 
must  be  excluded  from  the  response  during  [~T,  Oj  . We  normalized  to  unit  source 

strength  with  the  choice  P1°(t0)  = 1.  The  explicit  value  of  U+  is  fixed  by  the 

> 

choice  of  initial  conditions  r’(-r),  v(-t), 

Next,  a single  ion  is  permitted  to  evolve  from  t = 0 to  time  t,  when 

its  velocity  has  reached  the  asymptotic  value . No  further  interaction  with  the 

O 

external  electrons  is  included  hence,  we  choose  f%~0)  ~0^  and  P^1  (tQ  = 0) 

* l in  Eqn . (6a,  6b)  Conditional  probabilities  for  finding  the  atom  with  r(t) , 
v (t)  are  given  by: 


1/%,-r;  = Sx(o,t) 


15 


for  the  ionic  state  and 


-f 


v°(o,t)  = j d?-  Sjcoi'c')Q<t') 

■o 

tor  the  ground  state.  An  ion  may  reneutralize  at  any  time  after  t-  o.  Hence, 
V (O,  conto'?.  ‘-  integration  over  intermediate  times. 

The  single  atom  ejection  probabilities,!*' ° (-r. 0,1;,  for  the  ionic  and 
neutral  state,  follow  the  history  of  an  adsorbed  atom  which  responds  to  the 
external  beam  over  £~'"7'  ? 

Tc-r,  o,  t)  ■ U '('-t;  o)  y c o,  t) 

Ta(-r,0,t)  = UT(-T;0)  V°o,-C)  (8) 

These  quantities  are  closely  related  to  the  ionic  and  neutral  atom  desorption 

cross  sections  to  be  presented  below. 

Let  / denote  the  ionization  cross  section  for  an  adsorbed 

atom  at  the  point  r^  =T’’(t=0).  This  function  must  have  the  free  atomic  or 

molecular  ionization  cross  section  as  its  limit  when  z . Here.  £ 

o 

is  the  incident  electron  energy.  We  expect  ^ fa  ) to  be  linear  in  the 

upward  transition  rate  as  follows: 

)/r 


Here,  F is  the  incident  electron  flux. 


-r 


(9) 


This  result  is  recovered  after  first  noting  that  is  the  average 

ionization  rate  per  atom.  The  corresponding  average  ionization  cross  section  per 
atom  is  given  by: 


= -i—fc i-t' 

F-rJ 


rat:. » n . ; r v.  ^ , i . ...  c . *•*  * 
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In  the  limit  r~*0  + one  obtains  the  Instantaneous  response.  R ( t')  is 
smoothly  varying  for  ESD,  inasmuch  as  the  summations  over  final  state  electron 
coordinates  contained  in  R consume  any  delta  function  resonances.  One  factors 
R£(t'  = 0)  from  the  integration  above,  and  makes  the  replacement  S2  ( t',0)  ^ 

S2  ( -T|,0  ) s l to  obtain  Eq.  (9) 

The  cross  section  for  the  entire  process  resulting  in  ionic  desorption  of 
an  atom  Initially  in  the  ground  state  is  given  by: 


= K£(n)Vroti 


Here,  the  dependence  on  Rand  0 is  displayed.  The  step  function  was  inserted 
to  ensure  that  desorption  is  energetically  possible.  One  integrates  the  decay 
probability  appearing  in  the  exponent  from  point  to  point  as  the  trajectory  is 
followed.  The  asymptotic  velocity  generated  by  depends  on  the  choice 

°f  % ) &o'  One  chooses  t for  each  ion  path  so  that  z(t)  lies  far  from  the 


surface. 


The  analogous  ground  state  ejection  cross  section  is  given  by: 

<3 "'(ox)  = 

= ’XjnjEt'Qft'J  etpL- 

Ejri) -Eiff-’j+e/r')  - £,  <<*)) 

Here,  the  step  function  was  evaluated  at  the  reneutralization  polnt^ellminating 
trapped  particles.  E^  governs  the  path  integrations  during  | 0,t*y  , while  Ej 
determines  the  dynamics  during  j~t|  , t~J . Numerical  evaluation  of 


17 


is  expected  to  be  time  con  lining  due  to  the  integration  on  intermediate  times. 

Finally,  the  ensemble  averages  of  ejection  cross  sections  over  the 
phase  space  coordinates  is  iav.  - t*0.  The  formal  solutions  v/hich  yield  the 

total  ionic  and  neutral  atom  E .!>  t ; sections  are: 


(T  - \d  n d v0  ff' (h  OjXSh)  S/°,  t)e(^nX +£*(%)- 

J (iOa) 

^ t 

C r°  = j4,'%<d\  /f  °(r.,v.)  Xscr,)Cc \tt' S^o,t') Q (f) ■ 


■ 0(j:M\/0xt  - Et  (r'J  r £t  (*>)) 


(10b) 

The  equilibrium  distribution  Pj°  was  used  in  averaging,  rather  than  P^,  inasmuch 
as  the  ejection  cross  sections  are  already  linear  in  R.  Eq,  (10a)  for  (7”+  is 
employed  below  in  the  numerical  simulation  of  ESDIAD  for  0+  adsorbed  on 

Tungsten. 

Note  that  Eqn  (10a)  forCT+  is  a generalization  of  phenomenological 
expressions  which  have  appeared  in  the  literature.  Our  previous  work  on  one 
dimensional  Photodesorption  theory  produced  expressions  analogous  to  Eqns(lOa,  10b) 
for  the  optically  fluorescent  and  ground  state  cross  sections. 

The  result  expressed  in  Eqns  (10a),  (10b)  is  model  independent,  except 
for  general  restrictions  on  R.Q.E^,  and  P^°  outlined  above.  One  may  use 
Eqns{10  a,b)in  conjunction  with  model  surface  interactions  to  obtain  ESDIAD 
patterns  and  energy  spectra  within  the  Menzel-Gomer,  Redhead  ESD  picture. 

Finally,  note  that  the  transport  approach  is  applicable  to  a wide 
variety  of  surface  phenomena  involving  propagation  of  atcms.ions  or  molecules. 

In  particular,  the  cross  section  contributions  due  to  higher  order  excitation 
processes  in  ESD  may  be  obtained  without  difficulty. 


r 


IV  MODEL  CALCULATION  FOR  O*  DESORPTION  FROM  TUNGSTEN 


The  remainder  of  this  paper  applies  the  first  principles  kinetic  theory  in  a 
model  calculation  of  desorption  from  tungsten.  We  attempt  to  understand  the  interesting 
angular  and  energy  dependent  ESD  results  found  recently  on  W (111)6  by  numerically 
simulating  the  ESD  patterns  and  ion  energy  spectra  to  be  expected  for  atoms  adsorbed 
at  various  sites  on  the  surface.  Little  is  known  of  the  nature  of  possible  anti- 
bonding states  for  ions;  therefore,  simple  models  are  adopted  for  ion-solid  and 
atom-solid  interactions.  Absolute  cross  sections  were  not  calculated  due  to  lack 
of  information  about  wave  functions  on  surfaces.  We  concentrate  on  the  spectro- 
scopy of  ESD;  that  is,  in  this  paper  we  have  studied  only  the  polar  and  azimuthal 

angular  dependence  and  the  energy  spectra  of  emerging  ions. 

Experimentally,  a rich  diversity  of  ESDIAD  patterns  were  found  on 

W (111)  and  W (100).  For  high  coverage,  these  consist  of  groups  of  sharp  0+ 
emission  spots,  whose  geometry  is  a function  of  adsorption  and  annealing 
temperature.  0+  emission  normal  to  the  surface  was  found,  as  well  as  off 
normal  spots,  whose  number  and  orientation  clearly  reflect  the  underlying  sub- 
strate geometry. 

The  antibonding  state  potential  funnel  is  assumed  to  be  the  mechanism 

which  produces  focussed  spot  patterns.  For  a number  of  the  observed  ESDIADs, 

one  can  discount  the  possibility  that  off  axis  emission  is  due  to  thermally  acti- 

13 

vated  substrate  facetting.  Here,  only  flat  substrates  are  considered. 

The  peaks  in  ion  energy  spectra  were  found  to  occur  at  8 to  9 ev,  on 
6-8  5 

both  W(lll)  and  W(100).  This  agrees  with  earlier  studies  on  polycrystalline 
tungsten,  which  report  a spectrum  peaking  at  8.8  ev. 

The  qualitative  features  of  our  result  are  unaltered  by  the  following 
choice  of  models.  Spacial  variation  in  the  excitation  and  ^neutralization  rates 
was  neglected;  thus,  we  chose  arri  S2  to  be  constants.  When  vibrational 

states  are  highly  localized,  the  distortion  of  introduced  by  this  approxi- 

mation is  not  critical.  Ion  energy  spectra  will  be  skewed  toward  lower  energies, 
ar  i the  preferential  reneutralization  of  slower  ions  is  neglected. 


.<  L . ***  k 
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Near  a possible  adsorption  site,  we  modelled  E^  as  a parabolic 


potential: 


£t(\r) 


l 


Here,  "b  is  the  adsorption  site  coordinate.  The  chemisorption  energy  vo  does 
noc  appear  explicitly  m these  calculations.  The  spring  constant  tensor 


was  taken  to  be  isotropic,  for  simplicity.  One  fits 


to  the  best  experi- 


mental values. 


We  chose  for  P*  gaussian  distributions  both  in  position  and  velocity 
spaces..  This  choice  is  appropriate  for  particles  in  an  harmonic  potential 
well  both  for  high  temperature,  and  in  the  zero  temperature  case.  The  shape 
of  ion  energy  spectra  results  from  a folding  of  P°  against  in  these  cal- 

culations. However,  peak,  ion  energies  are  insensitive  to  the  width  of  P^  , 
The  sixfold  integration  on  phase  coordinates  r , vq  in  Eqn  (10a)  was 
accomplished  by  the  Monte  Carlo  method . 


In  modelling  E ions  were  approximated  as  point  charges  which 
respond  to  the  local  electrostatic  surface  potential.  The  character  of  the  domi- 
nent  anti -bonding  state  in  ESD  is  not  known;  for  simplicity,  we  chose  to  examine 
one  which  leaves  the  substrate  free  of  excitation.  is  taken  to  be  the  sum  of 
the  unexcited  substrate  Hartree  potential,  plus  an  image  contribution,  produced 
by  averaged  electron  density  fluctuations,  viz: 

Here,  q is  the  ionic  charge.  Bonding  effects  in  the  substrate  were  neglected. 

The  surface  Hartree  potential  was  represented  by  the  sum  of  free  tungsten 
atom  Hartree  potentials  centered  on  the  lattice  sites.  Similar  superpositions 

have  been  the  starting  point  for  self-consistent  band  structure  calculations  in 

. ..  * * 14-16 

bulk  tungsten. 
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17 

The  Herman-Skillman  wave  functions  were  used  to  evaluate  V 


numerically,  which  was  then  fitted  to  the  form: 

- /-29  r 

c*-)  = ?.  lx  e 


All  quantities  are  in  atomic  units.  As  a check,  we  evaluated  also  the  screened 

tungsten  atomic  potential  in  the  Fermi- Thomas  approximation.  At  realistic 

W-0+  separations,  the  resulting  potential  is  somewhat  smaller  than  the  above. 

18 ,19 

An  ion-surface  plasmon  hamiltonian  was  used  to  numerically  evaluate 
the  image  potential  Vj^  , for  a number  of  values  of  z and  r^,  the  plasma 
parameter.  The  result  for  singly  ionized  atoms  was  fitted,  to  the  following  form: 


- I 

+/f£  +/.3Vrs  r &Orsz) 


Image  charge  density  resides,  in  this  model,  on  an  approximate  electronic 
boundary  of  the  metal.  Including  plasmon  dispersion  shifts  the  divergence 
in  the  classical  result  at  the  dielectric  surface  to  the  interior  of  the  metal. 

Some  consequences  of  this  ESD  model,  on  W(lll)  are  now 
examined.  Extensive  numerical  computation  using  100°K.  for  the  desorption 
temperature  was  done,  at  a wide  variety  of  possible  adsorption  sites. 

The  most  striking  result  follows:  the  naive  model  constructed  for  the 
anti-handling  potential  strongly  focusses  ions  created  in  the  neighborhood  of  most 
plausible  candidates  for  adsorption  for  sites.  The  most  compact  simulated  O4 
emission  spots  are  associated  with  the  largest  mean  ion  energies.  Long  range 
Image  attraction  tends  to  defocus  the  spots.  However,  this  mechanism  is 
competitive  with  focussing  by  the  repulsive  component  of  E only  for  very  low 


energy  ions . 
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•~a.  h adr-  option  iu  candidate  produces  a single  emission  spot.  For  higher 

energy  ions,  the  local  value  of  y^-2-  at  l*ie  P°int  °*  ionization  provider 
the  domineut  smpuioo,  and  sols  ..nh  the  central  col  <titude  and  azimuthal 
orientation  eiative.  to  the  sub-  f ’ roi  the  resulting  emission  spot.  In  many 
cases  . the  do  ;u  ibutiou  to  L?  comes  from  i single  substrate  atom. 

Emission  -t  ‘ ots  the  chemisorption  bond  direction. 

V/e  were  successful  in  generating  simulated  ion  distributions  quite 
si'  liar  to  those  which  were  observed  on  W(lll).  Success,  however,  in  replica- 
ti  u the  visual  feature-  ’ ion  distributions  does  not  guarantee  that  this  model 
fci  the  antibonding  state  is  applicable.  In  particular,  the  simulated  ion  energy 
spectra  which  were  reduces,  peak  well  below  the  experimental  ion  energies. 

This  was  the  case  even  for  adsorption  site  candidates  whose  ESDIAD  pattern 
simulations  closeiy  resemble  those  observed.  This  finding  is  discussed  further 
in  the  concluding  section.  For  now,  we  remark  that  remodeling  of  the  anti- 
bonding state  is  required. 

Despite  the  limitations  of  this  model,  it  is  useful  to  explore  the 
most  plausible  picture  of  high  coverage  oxygen  chemisorption  it  can  generate, 
guided  by  experiment. 

Whenever  the  anti-bonding  state  potential  has  the  qualitative  behavior 
of  our  model,  the  following  rules  associate  the  gross  features  of  ESDIAD  spot 
patterns  with  adsorption  geometry.  Emission  spots  directed  off  the  surface 
normal  are  evidence  of  adsorption  at  a point  where  E^  is  assymmetric,  that  is, 
away  from  the  one  fold  degenerate  unit  cell  sites.  Here,  we  define  n fold 
degeneracy  to  mean  that  a particular  point  is  found  in  n equivalent  locations 
within  a single  unit  cell  whose  environments  are  identical  except  for  combinations 
of  rotation,  reflection,  and  translation.  For  example,  in  Fig.  (2)  all  the  points 
labelled  'a1 , '(5' , 1 ( ' , ' 1 are  degenerate.  At  the  high  symmetry  (one  fold 
degenerate)  sites,  the  anti-bonding  potential  possesses  approximate  azimuthal 
symmetry,  and  thus  leads  to  emission  peaked  normal  to  the  surface. 


Tl  j .i... . . »*  * 
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One  would  reconsider  this  interpretation  if  long  lived,  highly 
localized  substrate  excitations  could  be  shown  to  dominate  in  the  anti-bonding 
potential.  Such  excitations  might  be  localized  on  neighboring  adsorbate  or 
substrate  atoms,  and  could  result  in  non-normal  emission  of  atoms  adsorbed 
at  symmetric  points. 

Within  the  model  we  used,  symmetrical  groupings  of  off  axis  spots 
are  explained.  Assymmetric  points  in  the  unit  cell  are  several-fold  degenerate. 

An  adsorbate  domain  on  the  surface  in  which  such  sites  are  occupied  should 
be  found,  with  equal  occupation  probability,  in  all  degenerate  orientations. 

The  corresponding  ESDIAD  patterns  should  thus  contain  spot  superpositions 
suggestive  of  the  substrate  unit  cell  structure. 

We  illustrate  with  applications  of  the  model  to  W (111).  Fig.  (2)  shows 
the  W(lll)  surface  unit  cell.  Some  of  the  lengths  pertinent  to  the  W-O system 
in  are  tabulated  in  Table  1,  Here  we  are  guided  by  experiment.  Under 
conditions  for  which  facetting  of  the  surface  is  believed  not  to  occur, 
trios  of  0+  emission  spots  were  found,  having  the  substrate  orientation 
shown  in  Fig.  (2)  . These  were  found  in  conjunction  with  normal  emission, 
or  with  a second  spot  trio,  rotated  180°. 

In  this  simulation,  adsorption  at  the  highly  assymmetric  points  such 
as^1  or  'p1  in  fig.  (2)  can  produce  only  patterns  of  six  spots,  rotated  from  the 
observed  orientation.  Adsorption  directly  over  one  of  the  highly  symmetric 
points  'A' , ' B* , or  'C'  produces  only  emission  in  a narrow  cone  normal  to  the 
surface.  Adsorption  els.e where  along  the  unit  cell  mirror  planes  is  expected 
to  generate  spot  trios,  with  primary  oriented  spots,  as  in  Fig.  ^.corresponding 
to  adsorption  between  points  'A'  and  'C',  while  180°  rotated  trios  correspond 
to  occupied  sites  between  'C'  and  * B ' , or  'B'  and  'A'  . 


We  located  those  sites  which  yield  the  best  fit  with  this  anti-bonding 

potential,  to  the  experimental  ESDIAD  patterns  and  ion  energies.  A reasonable 

lower  limit  on  the  W-O  bond  length  was  chosen  as  l.95A°,  the  sum  of  Slaters 
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atomic  radii.  The  colatitude  of  . mulated  primary  lobes  was  fitted  to  the  ex- 
perimental^ lue  of  34°. 

The  most  plausible  site  candidate  found  for  the  origin  of  'primary' 
spot  trios  (i.e.,  those  with  the  orientation  shown  in  Fig. 2)  lies  almost  equi- 
distant from  the  'A'  and  'C'  layer  tungsten  atoms,  in  Fig.  (2),  along  the 
unit  cell  mirror  plane.  A site  almost  equidistant  from  the  'C‘  and  'B'  layer 
tungsten  atoms  is  this  models  best  choice  for  the  generator  of  180°  rotated 
spot  trios. 

For  both  of  these  sites,  the  mean  energy  of  emitted  ions  was  cal- 
culated to  lie  in  the  range  2 to  2.5  ev,  which  is  significantly  smaller  than  the 
experimental  value  of  8-9  ev.  For  the  normal  emission  spot  produced  by  ad- 
sorption directly  over  point  'C',  this  model  predicts  ions  energies  peaking  at 
1.3  ev.  Note  that  these  calculations  predict  peak  ion  energies  an  order  of 
magnitude  below  experiment,  for  adsorption  sites  directly  above  the  surface 
layer  tungsten  atoms  at  'A'.  Fig  (3)  shows  the  spot  patterns  for  sites  mentioned. 
V:  DISCUSSION,  CONCLUSIONS 

The  first  principles  kinetic  theory  presented  above  for  angular 
dependence  in  ESD  provides  a vehicle  for  testing  model  ion-solid  interactions. 
Calculations  performed  using  a simple  model  for  the  ground  state  surface  anti- 
bonding potential  predict  the  strong  focussing  of  ions,  like  that  observed  in 

+ + 

ESDIAD  experiments  involving  O emission.  Models  of  O chemisorption  on 
W were  inferred,  by  fitting  to  the  experimental  patterns.  Unfortunately,  poor 
agreement  between  our  calculated  ion  energy  distributions  and  the  experimental  ^ 
ones  weakens  the  plausibility  of  such  chemisorption  models,  for  high  oxygen 
coverage  on  tungsten.  However,  we  expect  these  model  calculations  to  apply 


to  low  oxygen  coverage  on  W for  which  is  believed  to  dissociate.  It 

is  likely  at  low  coverage,  that  the  metallic  substrate  will  support  only 

excitations  whose  lifetimes  are  much  shorter  than  ion  ejection  times. 

Complementary  model  calculations  for  desorption  from  W(100) 

11 

were  done.  Some  of  the  results  have  been  reported  in  the  literature.  Here,  we 
remark  that  the  ion  energy  distributions  obtained  in  model  computations  peak 

at  about  2.5  ev,  as  for  W(lll). 

Substantially  the  same  picture  of  adsorption  results  from  a range  of 

models  for  the  antibonding  potentials  (which  are  time  independent),  and  which 

Involve  basically  scale  and  range  changes.  More  careful  modelling  of  the 

potential  for  ground  state  surfaces  may  moderately  increase  the  calculated  ion 

energies.  For  example,  we  replaced  the  tungsten  wave  function  by  a set  which 
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includes  L-S  coupling,  then  repeated  calculations  for  some  selected  sites. 

Ion  energies  were  increased  by  a factor  of  roughly  1.5  to  2,  with  neglible  change 
in  the  simulated  ESDIA.D  patterns. 

It  appears  unlikely,  however,  that  a model  anti-bonding  state  in 
which  the  substrate  and  remaining  adsorbate  atoms  are  free  of  excitation  can 
adequately  explain  the  high  ion  energies  observed  for  pt  oxygen. 

It  is  interesting  to  note  that  ion  energy  spectra  peaking  at  1.7  ev, 
were  reported  for  ESD  of  H+  from  tungsten?  This  is  consistent  with  the  energy 
range  predicted  by  our  model. 

We  speculatively  suggest  that  excited  antibonding  states  involving 
the  surface  and  adsorbate  atoms  may  be  able  to  explain  the  large  0+  ejection 
energies  found  at  high  coverage.  Explicit  detailed  treatment  of  these  as  yet 
unknown  states  requires  cluster  calculations  beyond  the  scope  of  this  work. 

The  type  of  processes  we  suggest  are  analogous  to  those  involved 
in  electron  dissociation  of  free  molecules.  A collision  with  an  incident  electron 
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In  this  picture,  the  lifetime:?  of  surface  excitation  roust  be  long 
compared  to  ionic  flight  r<m«3S.  otherwise,  the  strong  mti-tonding  potentia 
shuts  off  immediately  after  its  creavicn.  leaving  ti  e ground  state  surface  po- 
tential to  provide  the  dominant  impulse  to  de sorbir.g  ions,  *s  n th;-  model  cal- 
culations presented  above,  I,ong  life  - Liras  for  sui f tea  excitations  may  be  justi- 
fied if  one  postulates  a transition  in  surface  behavior  from  metallic  to  insulating, 
as  coverage  approaches  a full  monolayer  for  oxygen,  "bote  that  the  experiment^  O 
energies  pertain  to  the  high  coverage  Pj  phase. 

Insight  into  the  physicat  reasons  fora  low  yield  P2  and  a high  yield 
p,  oxygen  phase  can  be  gained,  if  >hv  insulating  carrier  picture  is  a realistic 
one  for  high  coverage.  In  addition  to  supporting  .ong  lived  excitations,  such  a 
barrier  would  inhibit  ion  reneutralimtion.  Both  higher  ion  energies  and  smaller 
hopping  matrix  elements  would  appear  in  the  reneutralization  rate,  compared 
to  a low  coverage  surface.  Dramatic  increases  in  the  ion  survival  probability 
would  result,  due  to  the  exponential  dependence  of  S^on  Q(r). 

ESD  emission  of  H+  from  W can  probably  be  explained  using  the 
ground  state  surface  model,  with  no  need  to  postulate  the  formation  of  an  in- 
sulating barrier  at  high  coverage.  The  H+  yield  for  both  and  P2  hydrogen 
phases  is  small,  and  comparable  to  that  of  oxygen  in  the  low  coverage  phase. 
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We  can  devise  speculative  mechanisms  which  produce  geometrical 
spot  patterns  in  ESDIADs,  whether  it  is  the  ground  state  surface  potential  or  an 
excited  anti-bonding  state  which  governs  ion  emission.  If  the  anti-bonding 
state  is  associated  with  highly  localized  excitations,  off  axis  spots  can  result 
even  in  the  case  where  adsorption  occurs  at  the  highly  symmetric  points  of  the 
unit  cell.  Note,  for  example,  on  Fig.  (2),  that  an  adsorbate  ion  created  above 
point  'C',  together  with  excitation  of  an  adjacent  'A'  layer  tungsten  atom,  can 
proauce  downward  facing  spot  trios.  If  excited  states  are  multi-centered,  the 
anti-bonding  potential  will  have  essentially  the  same  symmetry  as  the  ground 
state  potential.  Off  axis  spot  groups  then  correspond  to  adsorption  at  points 
in  the  unit  cell  away  from  the  symmetry  sites. 

The  speculative  nature  of  some  of  our  remarks  indicates  the  need 
for  accelerated  theoretical  study  of  ESD  mechanisms,  after  which  unambiguous 
adsorption  site  determinations  using  ESDIAD  may  be  made. 
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FIGURE  CAPTIONS 

Fig.  1 The  one  dimensional  qualitative  behavior  of  effective  potentials  for 
ionic  ESD,  as  functions  of  z coordinate  . V is  the  chemisorption 
energy.  The  z coordinate  of  the  adsorption  site  is  'b*.  The 
The  metal  + adsorbate  system  in  its  ground  state,  designated  by 
(M  + A),  has  the  potential  function  EjGr).  Ai|  anti-bonding  potential 

E (r)  leading  to  ionic  ESD  is  denoted  (M  + A ).  The  ^u.b^trate 
2 * | 
here  Is  free  of  excitation.  The  dashed  line  (M  + A ) suggests  the 
behavior  of  the  effective  potential  for  an  ion  interacting  with  the 
excited  substrate.  The  asymptotic  energy  difference  E - Ej  is 
V.,  the  atomic  ionization  potential,  if  both  incident  ana  ejected 
electrons  are  Auger  emitted  to  the  vacuum. 

Fig.  2 The  W(lll)  surface  unit  cell.  Large  spheres  depict  tungsten  atoms. 

The  'A*  atoms  lie  in  the  surface  layer.  'B*  and  ‘C’  atoms  reside  in 
the  second  and  third  planes  in  the  solid.  The  A- A distance  is 
8.44  a . Dashed  lines  through  “B',  ’C‘  are  bilateral  symmetry 
axes.  °' Primary'  ESD  emission  lobes  for  O on  W(lll)  have  been 
observed,  with  their  substrate  registry  as  indicated  by  the  trio  of 
heavy  circles. 

Fig.  3 Simulated  ESD  spot  patterns  for  low  coverage  of  Qxycjen  on  Wftll).  Fig. 
3a  shows  the  ion  distributions  for  adsorption  1.95  A . over  point 
'C'  on  Fig.  2.  Fig.  3b  corresponds  to  adsorption  almost  equidistant 
from  tungsten  atoms  at  'A'  and  'C',  while  Fig.  3c  pertains  to 
adsorption  between  sites  ‘C‘  and  'B'.  This  series  provides  the 
best  fit  to  experimental  patterns  of  Madey,  et  al. 
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